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A description of spin Faraday rotation, Kerr rotation and ellipticity signals for single- and multi- 
layer ensembles of singly charged quantum dots (QDs) is developed. The microscopic theory consid- 
ers both the single pump-pulse excitation and the effect of a train of such pulses, which in the case of 
long resident-electron spin coherence time leads to a stationary distribution of the electron spin po- 
larization. The calculations performed for single-color and two-color pump-probe setups show that 
the three experimental techniques: Faraday rotation, Kerr rotation and ellipticity measurements 
provide complementary information about an inhomogeneous ensemble of QDs. The microscopic 
theory developed for a three-dimensional ensemble of QDs is shown to agree with the phenomenolog- 
ical description of these effects. The typical time-dependent traces of pump-probe Faraday rotation, 
Kerr rotation and ellipticity signals are calculated for various experimental conditions. 

PACS numbers: 78.67.Hc,78.47.-p,71.35.-y 
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I. INTRODUCTION 

It is impossible to overestimate the role, which pump- 
probe spin Faraday and Kerr rotation measurements have 
played and continue to play in developing of spintronics, 
a new area of science and technology that tries to utilize 
an electron spin, in addition to its charge, in various semi- 
conductor devicesii^. The discovery of a very long spin 
coherence time in bulk GaAs^ and II- VI compound quan- 
tum wells^ using these highly sensitive techniques was 
one of the cornerstones for the initiation of spintronics, 
and today pump-probe spin-dependent spectroscopy has 
become a common way to study carrier spin coherence in 
bulk crystals^, quantum wells (QWs)ii^ and quantum 
dot (QD) samples^iiiSiiiii^ii^. At the same time Kerr ro- 
tation measurements have become the most visual and 
impressive method to study electron spin transpor t^^i^^ , 
spin accumulation and injectio n^^i^^ , and the spin-Hall 
effecO^^. 

The schematic illustration of the pump- probe measure- 
ment techniques is shown in Fig. 1. A first short intense 
pulse of circularly polarized light (a pump pulse) gener- 
ates the nonequilibrium spin-oriented electrons and holes 
and creates a macroscopic spin polarization^^. In a con- 
stant transverse magnetic field, B , applied to the sam- 
ple the macroscopic polarization starts to precess around 
the field direction. On the microscopic single-spin level, 
the precession is connected with a coherent superposi- 
tion of two spin levels split by the magnetic field. The 
superposition is created by a short pulse of circularly po- 
larized light and the quantum mechanical beating of this 
coherent superposition occurs at the Larmor precession 
frequency of the applied magnetic field, = gefJ-BB/h , 



where ge is the electron g -factor, and hb is the Bohr 
magneton. The optically created polarization and its pre- 
cession can be probed by short pulses of linearly polar- 
ized light via rotation of their polarization plane after 
the propagation through the photoexited medium (Fara- 
day effect) or reflection from this medium (Kerr effect). 
The short probe pulses of linearly polarized light show a 
remarkable sensitivity to the practically instant popula- 
tion of electron and hole spin sublevels. The pump-probe 
techniques allow one to study spin dynamics of resident 
carriers that are also polarized by the pump pulse during 
their coherence time, which exceeds the typical photolu- 
minescence decay time by several orders of magnitude. 
These advantages make the pump-probe Faraday and 
Kerr rotation techniques suggested more than 15 years 
g^gQ 2i^22^23,24 -j-q j^^^ powerful tool to study the carrier 
spin dynamics. 

The pump-probe spin-dependent rotation techniques 
are especially useful for manipulation and measurement 
on electron spin polarization in singly charged QDs due 
to a very long coherence time of resident-electron spins. 
This time in QDs could be as long as several microsec- 
ondsiii^^ and it exceeds the spin coherence time mea- 
sured in bulk GaAs^i2^ by two orders of magnitude. The 
resonant short pulse excitation of such QDs by circu- 
larly polarized light leads to practically deterministic cre- 
ation of an electron spin polarizatioi>i^. A train of such 
pump pulses results in complete synchronization of elec- 
tron spins if their precession frequency satisfies the phase 
synchronization condition (PSC), which is fulfilled when 
the train repetition period being shorter than the sin- 
gle electron spin coherence time is equal to an integer 
number of the Larmor precession periodsSI. This syn- 
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FIG. 1: Schematic illustration of (a) the pump-probe Faraday 
and Kerr rotation measurement techniques and (b) a coherent 
superposition created by a short pulse of cr+ polarized light 
from the two spin states split in a transverse magnetic field, 
B . Of and 6k are the Faraday and Kerr rotation angles, 
respectively. 



chronization leads to the mode locking of electron spin 
coherence in an ensemble of QDsii, the effect which al- 
lows one to overcome the dephasing of electron spin po- 
larization connected with a dispersion of precession fre- 
quencies and to control the ensemble polarization during 
a single electron spin coherence time. The mode-locking 
effect in self-organized (In,Ga)As/GaAs QDs is enhanced 
by the nuclear induced frequency focusing effect, which 
shifts the electron spin precession frequencies of the en- 
tire ensemble of QDs to the modes satisfying the PSC— . 
This last phenomenon opens exciting opportunities to 
create an entire ensemble of QDs with more than one 
million spins having a single precession frequcncjs2^ and 
their controlled manipulation by short pulses of various 
polarizations^Si^. 

Despite the tremendous success of the pump-probe 
Faraday and Kerr rotation measurement techniques, 
their microscopic theoretical descriptions for QD struc- 
tures are absent to the best of our knowledge. For 
bulk semiconductors and QWs the theory of magneto- 
optical photoinduced Faraday effect was developed, re- 
spectively, by Aronov and Ivchenko^^ and by the au- 
thors of Rcfsiii^^i^. In the present paper we develop 
such a theory for the array of QDs considered as an en- 
semble of independent localized oscillating dipoles. This 
approximation, which generally imposes a restriction on 
the QD concentration, is usually sufficiently accurate in 
self-organized QD samples like those used, e.g., in RefJ^. 

In the standard pump-probe Faraday rotation experi- 
ments, the transmitted probe light is split into two lin- 
early polarized beams with orthogonal polarizations ori- 
ented at ±45° angles relative to the initial light polar- 
ization (or the polarization of the probe light transmit- 
ted through the unpumped sample). Then the difference 
of the time-integrated intensities of the split beams is 
measured as a function of the delay between pump and 
probe pulses^. To describe the experimental setup we 
introduce two pairs of axes, x,y and x',y' , rotated by 



a 45° angle with respect to each other. The initial po- 
larization of the probe light is along x axis. This allows 
us to define the experimentally measured spin Faraday 
signal, , as 



J- — lim 

Z — ' + 00 



dt, (1) 



where E^,'{z,t) and Ey}{z,t) are respectively the x' - 
and y' -components of an electric field of the transmit- 
ted probe light at time t . They are connected with 
the X- and y -components of the electric field: E^' = 
{E^ - Ey)/V2 and Ey, = {E^ + Ey)/V2. Equation 
([1]) assumes the probe light source to be positioned at 
z — > — oo . 

In this form Eq. ([1]) is derived for the case where the 
sample is subject to a periodic train of pump and probe 
pulses repeated with a certain period Tr . The integra- 
tion in Eq. ([T]) takes place over the measurement time, 
Texp , which exceeds by far all other time constants in 
the experiment, such as spin precession and relaxation 
times and pulse repetition period. For the case of a sin- 
gle pump and single probe pulse the integration in Eq. ^ 
is effectively carried out during the probe pulse duration, 

The Kerr effect is measured in the reflection geometry 
and its magnitude is defined as 



lim 



where E^,\z,t) and E^y,\z,t) are respectively the x' - 
and y' -components of the reflected probe pulse electric 
field. The probe-pulse rotation angles in the Faraday 
and Kerr rotation measurements can approximately be 
expressed in a simple form^S, 



|i?«(.,t)p-|i?W(,,t)|^ di, (2) 



a:/2x. 



(3) 



if I^f.kI <C 1 . Here I is the total time- integrated inten- 
sity of the transmitted or reflected probe pulse, respec- 
tively. 

The pump-probe ellipticity measurement is another 
way to study spin dynamics in samples photoexcited by 
the polarized pump pulso^i. The experimental setup in 
this case is similar to setup used for the Faraday rotation 
measurements but with 1/4 waveplate. The ellipticity 
signal in transmission is deflned by the following inte- 
grated difference 



E 



lim 

; — ^+oo 



\Ef^{z,t)r dt, (4) 



where e'^^1 = {E'i^ T \E^y^)/^/2 are the circular cr+ 
and components of the transmitted probe pulse. For 
small ellipticity values, the so-called angle of ellipticity, 
£ , is given by 



£/2I. 



(5) 
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In this paper we calculate microscopically the magni- 
tudes of the single- and two-color pump-probe Faraday 
and Kerr rotation signals, T and K, , and the elliptic- 
ity signal £ in an ensemble of singly charged QDs reso- 
nantly excited by a single short light pulse of an arbitrary 
shape or an infinite train of such pulses. The calcula- 
tions show that the three measurement techniques ex- 
plore spin polarization properties of different subsets of 
QDs due to inhomogeneous broadening of the resonant 
transition energies in the QD ensemble, and the results 
of these measurements may be nonidentical. As a result, 
the time-dependent traces measured by these three tech- 
niques significantly differ from each other and are very 
sensitive functions of the pump and probe excitation fre- 
quencies, their detuning, and the dependence of the elec- 
tron g -factor and the oscillator transition strength on 
the trion excitation frequencies. The traces of the Kerr 
rotation signal depend also on the thickness of the cap 
layer. The average electron spin precession frequencies 
measured in these experiments can differ as well. 

Our paper is organized as follows. In Sec. [Tllwe pro- 
vide a theoretical description of electron spin polariza- 
tion created in an ensemble of QDs by a single short 
pulse of circularly polarized light or by an infinite train of 
such pulses. The general microscopic theory of probing 
this spin polarization in the QD ensemble is presented 
in Sec. IIIII For a three-dimensional ensemble of QDs 
we compare the developed approach with the standard 
phenomenological description of magneto-optical pump- 
probe effects within the effective-medium approximation. 
In the final Sec. IV we calculate typical time-dependent 
traces of the two color pump-probe Faraday and Kerr ro- 
tation and the ellipticity signals and briefly compare the 
obtained results with available experimental data. 



II. CREATION OF ELECTRON SPIN 
POLARIZATION IN SINGLY CHARGED 
QUANTUM DOTS 

In what follows we consider a planar array of singly 
charged zinc blende based QDs grown along the axis 
z 11 [001] . The QDs are self-organized during molecu- 
lar beam epitaxy growth on the wetting layer. The lat- 
eral size remarkably exceeds their height, which serves a 
quantization axis for the electron ±1/2 spin states and 
heavy- hole ±3/2 spin states responsible for the domi- 
nating optical transitions. In the absence of a magnetic 
field, the ground state of a singly charged QD is two-fold 
spin degenerate. The first excited state of such a QD un- 
der interband transitions is a singlet trion, which consists 
of two electrons occupying the same size-quantized level 
with opposite spins and a heavy hole in one of the two 
degenerate states: ±3/2 . The optical selection rules for 
the resonant excitation of these trion states and their ra- 
diative decay are very restrictive. The ±3/2 trion states 
can be created only by the cr"*" circularly polarized light 
and only in a QD where the resident electron has the 



spin projection ±1/2. These ±3/2 trion states can ra- 
diatively decay only into the initial ±1/2 spin states. 
At the same time (7+ circularly polarized light does not 
affect an electron with the spin projection —1/2. The 
same rules with the sign reversal "±" " — " are ap- 
plied for the optical excitation of the —3/2 trion state. 
It is important to notice that the singlet trion in these 
QDs does not have an optical transition dipole compo- 
nent along the z axis. 

A. Pumping of electron spins in quantum dots 

Let us first consider the effect of QD photoexcitation 
by a short electromagnetic pulse with the carrier fre- 
quency Wp close to the trion resonant frequency wq ■ 
We also assume that the pulse duration time Tp is short 
compared with other times: the spin relaxation times of a 
resident electron and a photohole forming the trion; the 
trion radiative lifetime; and the spin precession period 
of an electron and a heavy hole in an external magnetic 
field. According to the selection rules the interaction of 
the QD with the electromagnetic wave is described by 
the Hamiltonian 

I>(t) = - y [rf+(r)£;,+ (r, t) ± d_ [r)E,- (r , t)\d?r , (6) 

where d±(r) — Xd^iy) ± i(ij,(r)]/-\/2 are the circularly 
polarized components of the dipole moment density op- 
erator, and E^±{r,t) are the circularly polarized com- 
ponents of the electric field of a quasi monochromatic 
electromagnetic wave. The electric field of this wave is 
defined as 

E{r, t) = E^+ (r, t)o+ + E„- (r, t)o_ ± c.c. , (7) 

where o± are the circularly polarized unit vectors re- 
lated to the unit vectors Ox \\ x and Oy \\ y by 
o± = {ox ± iOy)/y/2 . Here the both components E^+ 
and E^- are proportional to the exponential function 

The incident electromagnetic field induces optical tran- 
sitions between the electron state and the trion state cre- 
ating a coherent superposition of them. In accordance 
with the selection rules the circularly polarized light 
creates a superposition of the ±1/2 electron and ±3/2 
trion states while the a~ polarized light creates a super- 
position of the — 1/2 electron and —3/2 trion states. In 
order to describe these superpositions it is convenient to 
introduce a four component wavefunction 

*= (^'1/2, -0-1/2, ■03/2, ■0-3/2) , (8) 

where the ±1/2 subscripts denote the electron spin pro- 
jection and ±3/2 refer to the spin projection of a hole 
in the trion. The electron spin polarization is expressed 
in terms of 0±i/2 as follows 

S. = (|V'l/2|'-|V'-l/2|')/2, 

Sx = Re(0i/2V'* 1/2) , = -Ini(0l/20* 1/2) ■ (9) 
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Hereafter we completely neglect all other excited states of 
a QD, e.g., triplet trion states, and treat the QD optical 
excitation within the four-level model. In this approx- 
imation the action of a short pulse on the charged QD 
can be described by the following equations 



y+(t)Vl/2, (10) 

1 

hiJoij-3/2+V-{t)i:_l/2 , (11) 



ihlpi/2 = (0^3/2 
ihlp_3/2 
ihlp_i/2 = Vl{t)l(;_3/2 



Here -0 = dij)/dt and the time-dependent matrix ele- 
ments V±{t) = — j 6{r)E^±{r^t)<l^r describe the light 
interaction with a QD. The strength of this interaction 
is characterized by the effective transition dipolo^ 



d(r) = (l/2|d_(r)|3/2) = (-l/2|d+(r)|-3/2) 



Wo mo 



-F(r,r) 



(12) 



which is the matrix element of the operators d±{r) in 
Eq. IH]) calculated between the wave functions of the va- 
lence band, I ± 3/2) , and the conduction band, |±l/2), 
all taken in the electron representation. In Eq. p^ . 
e is the electron charge, mo is the free electron mass, 
and pcv = (S|px|X) = (S|p.y|Y) = (S|p^|Z) is the inter- 
band matrix element of the momentum operator taken 
between the conduction- and valence-band Bloch func- 
tions at the r point of the Brillouin zone, S and (X , 
Y , Z) , respectively. Finally.— 

F{r,,ru) = Mrh)v^e'\r,) U\' ^,{r')^^^'\r'), (13) 



where (f^e'^^ and i^h are, respectively, the electron and 
heavy-hole envelope functions in a trion, and ip^, is the 
envelope function of a single (resident) electron confined 
in a QD. These wave functions are chosen to be real. 
In derivation of Eq. we assumed that the QD is 

small enough so that the electron and hole motion in a 
trion can be treated independently, but allowed for the 
different orbital functions of the resident electron and 
the electrons in a singlet trion4S Note, that in Eq. 
F(re,rft,) is taken at the coinciding coordinates of the 
carriers, re ~ Vh ~ r . 

Due to a very short timescale of the pump pulse we 
completely neglect the electron and hole spin precession 
during pulse action as well as spin dephasing and radia- 
tive decay processes in Eqs. pU)) and pTjl . The dynamics 
of the electron spin polarization described by these equa- 
tions was considered before in Ref.— in the case of a rect- 
angular shape pulse with the resonant carrying frequency 
tjp = ujQ . This paper generalizes the consideration for 
detuned pulses of an arbitrary shape. 

To be specific we consider the excitation of the QD 
by a polarized light pulse; the difference in the dy- 
namics of the electron spin polarization created under the 



a~ photoexcitation is briefly discussed below. Before the 
pump pulse arrival the QD is always in the ground state 
because the pump repetition period is much longer than 
the trion lifetime in the QD. For the polarized exci- 
tation the component "0-1/2 of the QD wavefunction is 



conserved, and -0- 
set of Eqs. 



3/2 



. This allows one to reduce the 



To)) to a single equation for the component 







1/2 




^'1/2 + /'(i)^l/2 = . 



(14) 



Here lu' = lUj, — ojq is the detuning between the pump 
frequency and the trion resonance frequency, and f{t) is 
a smooth envelope of pump pulse defined as 



fit) 



dir)E,^ (r,i)d3 



It follows from Eq. p4p that the values 0i/2 and 
"01/2 before and after the pulse action are connected 
linearly^^ . Taking into account that the initial conditions 
for Eqs. pU)) are ■0i/2(— oo) = const, 03/2(— oo) = 
and, therefore, 0i/2(~oo) = , one can represent the so- 
lution of Eq. p4|) at t ^ Tp , i.e., after the pulse is over, 
as 



V'i/2(oo) = (3e'*0i/2(-oo) . 



(15) 



Here the real coefficient Q satisfies the condition ^ 
Q ^ 1 and the phase <& can be chosen in the interval 
between — tt and tt . Both parameters are determined 
by the pump pulse shape, power and detuning. Equa- 
tions ([9]) and p5|) determine the modification of the elec- 
tron spin polarization from a short pulse of an arbitrary 
shape. Taking into account that ip-1/2 is conserved un- 
der the polarized excitation, the electron spin before 
the pulse arrival, = (S^ , Sy ,S~) , and just after the 
end of the pulse, ~ (5^, Sy,S^) , are connected by 



St 

Sx 

s+ 



1 



-S: 



4 2 

Q cos $5"" + Q sin $5"^" 
Q cos ^Sy - Q sin $5"^ 



(16a) 

(16b) 
(16c) 



Although Eqs. PB)) are derived for the pure spin states 
they are valid as well for the mixed states whenever the 
pulse duration Tp is much shorter than the spin relax- 
ation times in the QD and the time of electron and hole 
spin precession in a transverse magnetic field. 

Using Eqs. pU]) one can show that |'0i/2(i)P + 
lV'3/2(i)P = |V'i/2(— oo)p . It follows then that the z 
component of the post-pulse trion spin polarization for- 



mally defined as Jz 
equal to 



(l^3/2(«3)p-|0_3/2(oo)P)/2 



IS 



Jz — 



S ^ 



(17) 
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Derivation of the relation between S'+ and S~ estab- 
hshed by the cr~ circularly-polarized pulse gives equa- 
tions similar to Eqs. (fT6|) . One has, however, to change 
the sign of the first term in Eq. (|16ap and replace $ by 
-$ in Eqs. (|16bp and (fTec]) . 

Equations (fT6|) and p7|) are the main result of this 
Section. They show how the short pulses of circularly 
polarized light create and control the electron spin po- 
larization in n -type QDs under the resonant trion ex- 
citation. At low pump intensities the electron spin is 
weakly affected and a value of Qe'* slightly deviates 
from unity. For high intensities the coefficient Q notice- 
ably decreases and the phase $ shifts from zero. One 
can see from Eq. ()16ap that the cr+ circularly polarized 
pulse modifies the z -component of electron spin polar- 
ization by S+ - S- = (Q2 - 1)(1 2S-)/A . The pump 
pulse also leads to the in-plane rotation of the electron 
spin polarization along the light propagation direction for 
pulses with $ 7^ (see Eqs. (jl6bp and (|16cp ) similarly 
to spin rotation in a longitudinal magnetic field B |j z . 

It is easy to show that the phase $ becomes nonzero 
due to detuning between the resonant and pump frequen- 
cies. Indeed, under resonant conditions, to' — , or small 
detuning, wVp <^ 1 , one can neglect iuj' as compared 
with f{t)/f{t) ^ in Eq. (fTH) and produce the solu- 
tion in the form 



V'l/2(0 = l/'l/2(-00)c0S 



fit'W 



(18) 



The direct comparison with Eq. pS]) gives <& = and 
Q = cos (6/2) , where 



e 



/(t')dt' 



(19) 



is the effective pulse area. This consideration shows that 
only detuned pump pulses give rise to $ ^ and gen- 
erate an effective magnetic field acting on an electron 
spin in the QD. The pulse tuned resonantly to the trion 
transition causes no rotation of the in-plane spin compo- 
nents, and the electron spin dynamics is independent of 
the pulse shape and is controlled only by the pulse area 

e. 

Let us now analyze the dependence of Q and $ on 
the pulse parameters for the detuned pulses. In the gen- 
eral case of an arbitrary pump pulse power, its arbitrary 
detuning and shape, Eq. (jl4p can be solved only numer- 
ically. The analytical solutions of Eq. p4|) are known 
for the two pulse shapes: (i) for pulses with rectangu- 
lar shape, f{t) = /o = const 7^ for \t\ < Tp/2 and 
f{t) = otherwise, (ii) for smooth pulses of the shape 
suggested by Rosen and Zener— : 



fit) = 



cosh (TTt/Tp) ' 



(20) 



where the coefficient /i is a measure of the pulse electric 
field strength. The effective areas of these pulses are 
equal to 9 = 2/orp and 9 = 2yLtTp , respectively. 



For f{t) taken in the form of Eq. ([20)) one can write 
the solution of Eq. fH]) following Ref4^ in terms of the 
hypergeometric function 



1pl/2{t) = -01/2(-Oo) X 

"9_ _ 9_ 1 

2^'~2^' 2 ' 



2F1 



■ tanh 





1' 


(?) 


+ 2 



(21) 



where 9 = 2/irp and the dimensionless detuning y = 
lu'tp/(2tt) . This leads to explicit expressions for Q and 



Q 



iy) 



sin^(9/2) 
cosh^ (Try) 



arg- 



r(i-f -iy)r(i + £-i,) 

In the case of a rectangular shaped pulse we obtain 



(22) 
(23) 



Q 



92 

X 



— sm - 



ny - (I), 



(24) 
(25) 

(26) 



where 9 = 2/oTp , effective Rabi frequency 



and sin = (y/Qx) sin(x/2) . One can see from Eqs. 
and ((25)) that $ changes its sign with reversal of the 
detuning parameter y . For circularly polarized pump 
pulses, this means that the sign of detuning determines 
the rotation direction of the electron spin polarization 
and its reversal is similar to switching the effective mag- 
netic field from one direction to the opposite. 

According to Eqs. (|16p and (fT7|) the reorientation of 
electron spin polarization by short pulses with Q ~ 
or Q = 1 becomes deterministic and can be used for 
controllable manipulation of a single electron spin. The 
pump pulses with Q = completely erase, for any de- 
tuning, the in-plane spin component a^'^i'^^ leading to the 
alignment of the electron spin along the z axis. On the 
other hand, the pulses with Q = 1 lead to a control- 
lable rotation of the electron spin in the {x, y) plane 
and make no effect on the spin z component. The ro- 
tation angle equals to $ and is determined by the de- 
tuning^l. It follows from Eqs. (I21])-(I1S1) that the value 
Q = can be reached only with the pulses tuned to the 
resonance ( y = ) and having the areas 9 = tt, Stt, ... 
(the so-called tt -pulses). In the case of Rosen&Zener 
pulses^ the condition Q = 1 is realized for any detun- 
ing if 9 = 27r,47r, ... ( 27r -pulses). For the rectangular 
pulses this condition can be reached only for certain com- 
binations of the detuning and pulse area when the effec- 
tive Rabi frequency x = ^ (27r?/)2 + 9^ = 27riV with N 
being an integer. 
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Figure[2]shows the calculated dependences of Q and <!> 
on the detuning for four pulse areas , each for rectan- 
gular and Rosen&Zener pulse shapes. For large detuning, 
y = \lo^—ujo\tp/2tt 1 , Q is close to 1 , $ tends to . 
Therefore, the electron spin state in a QD is unaffected 
by the strongly detuned pulses. For the area 9 = tt , the 
function Q{y) has a sharp dip at y = reaching zero 
value at this particular point. Thus, for tt -pulses tuned 
to the resonance ut^ the parameter Q vanishes and, as 
stated above, the pump pulses suppress the transverse 
spin components and Sy . Deviation of O from tt 
converts the dip into a smooth minimum. In accordance 



S~ = ; before the pulse arrival. On can see that the 
rectangular pulses with 8 > tt result in intensive oscilla- 
tions of 5+ at large-scale detuning, whereas oscillations 
are completely absent for the Rosen&Zener pulses. The 
oscillations are again connected with the shape of the 
Fourier transform of a rectangular signal. One can see in 
Fig. E^d) that detuned pulses of a different shape create 
a completely different electron spin polarization. 



B. Temporal dynamics of electron and trion spin 
polarization after the pulse 



Rosen & Zener 
pulse 
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FIG. 2: Dependence of Q and $ on detuning y = (cjp — 
LJo)T"p/27r for pulses having the rectangular shape (a,b) and 
Rosen&Zener shape (c,d) calculated for several pulse areas 
Q = 7r/2, TT, 37r/2, I-k . Insets demonstrate shapes of pulses. 



with Eq. (im, the pulses with the areas O and 27r — O 
produce the same Q . The signs of the phase <I> and 
detuning are opposite, and makes, at zero detuning, 
an abrupt jump from its positive maximum value to the 
negative minimum value. The altitude of this jump rises 
along with the pump pulse area. Although the both pulse 
shapes lead to generally similar dependences of Q or <I> 
on the detuning, the rectangular shape pulses give rise to 
additional oscillations on the curves Q{y),^{y) clearly 
seen in Figs. [^Ja), [^JJb). These oscillations are connected 
with the oscillating character of the Fourier transform 
of a rectangular shaped signal. They are absent for the 
Rosen-Zencr pulse, for the general case of a smooth pulse 
these oscillations are much weaker than for the rectan- 
gular pulse. One can see that, for 6 7^ tt , Q{y) has 
several minima at y ^ Q . This occurs because at these 
detuning the effective Rabi frequency, x , approaches to 
the Ntt with N = 1,3,5... . 

Figure [3] compares the z -component of electron 
spin polarization created by rectangular (dash-dot black 
curves) and Rosen&Zener (solid red curves) pulses of dif- 
ferent areas O in singly charged QDs with zero spin. 



The temporal dynamics of electron and trion spin po- 
larization in a QD after its excitation by the short pump 
pulse can be determined from the kinetic equations for 
the electron spin polarization S and z -component 
of the trion spin polarization— i^'^i^^ 



. J. 



TQD 



s 

A 

Ts.h 



TQD 



(27) 



The equations take into account (i) the precession of elec- 
tron spins in the in-plane magnetic field B with fre- 
quency rj , (ii) the electron spin relaxation, (iii) the spin 
relaxation of a hole in the trion, (iv) the radiative decay 
of the trion, and, finally, (v) the partial suppression of 
the electron spin polarization created by the pulse after 
the trion recombination. It is ignoring, however, the hole 
precession. Here r^^e , TsJi and tqd are the single elec- 
tron spin relaxation time, the hole spin relaxation time 
and the trion lifetime, respectively. 

The detailed analysis of the spin dynamics in the 
coupled electron-trion system has been carried out in 
Refsjii^l (see also^i^). It is instructive to consider the 
simplest case of these dynamics when the hole-in-trion 
spin relaxation is much faster than the trion radiative 
lifetime, Ts.h <C tqu , or when the electron spin preces- 
sion time 27r/il <C tq£i . In the both limits the elec- 
tron remaining in the QD after the trion recombination 
becomes completely depolarized and its contribution to 
the electron spin polarization created during the pulse is 
completely negligible^iiii^. In this case the precession of 
electron spin polarization in a transverse magnetic field 
after the trion decay is described by the following set of 
equation! 
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Sz{t) 

Sy{t) 



[S'+cosm + S'+sinm]e" 
= [5+ cos nt - 5"+ sin nt]e 



(28) 



where time t is referred to the end of excitation pulse, 
and the electron spin polarization components {a = 
x,y, z) created by the pulse are defined by Eqs. ([TB)) . 
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FIG. 3; Effect of detuning u-p—ujo on tlie z; -component of electron spin poiarization, , created by tfie rectangular (dasli-dot 
black curves) and Rosen&Zener (solid red curves) pulses. Curves are calculated using Eq. (|16ap with 5*^ = and Eqs. (|22l 
I24[) for the pulse areas: Q = O.Stt (a), tt (b), I.Stt (c), and 2-r (d). 



Electron spin polarization created in QDs by an 
infinite train of short pulses 



In the pump-probe Faraday and Kerr rotation experi- 
ments the sample is usually subjected to a train of pump 
pulses that follow with a certain repetition period Tr . 
If the time is comparable with or smaller than the 
single electron spin relaxation time in a QD, < Tg^e , 
the electrons retain the memory of being exposed to the 
previous pulses. The infinite train of pulses creates a 
steady state of the electron spin polarization in the QDs 
periodically varying in time with the same period , 
which leads to a number of remarkable phenomena such 
as the resonant spin amplification^!^ and the mode lock- 
ing of electron spin cohcrenceiiii^. The time evolution 
between the pulses is described by Eq. ([28]) where 
(a = x,y,z) are the components of spin polarization 
taken at the end of any pump pulse. To find these com- 
ponents one should associate the polarization at the 
moment t = Tr with the spin polarization at the mo- 
ment before arrival of the next pulse, S{Tjj) = , and 
interconnect it with according to Eqs. pG)) . As a 
result we obtain self-consistent equations for the compo- 
nents . Solving them and substituting the solution 
into Eqs. (|16p we find the components which can be 
written as 













4A 




1-g 




4A 



-e sin(f]TR) , 



-e X 



Q(cos<i? - A'sin$)e' 



(29) 



where 
A 



Q(l + Q^) .-2TWr. 



Q(cos $ - K sin $) cos(r2TH) 
(cos $ - A' sin $) , 



K 



1 - Qe-"^"/^-." COS* 



(30) 



One can check that, for a periodic train of pulses of 
arbitrary intensity and shape, the electron spin polar- 
ization reaches the highest value at the magnetic field 
satisfying the PSC condition = 2ttN/Tr , where N 
is an integeri^iiiii^d^. For such electrons cos {^Tu) = 1 
and the z -component of their spin polarization at the 
moment of pulse arrival can be written as 



1 - 



(31) 



Equation (|3ip shows that the maximum value of \S^\ is 
independent of the phase shift $ between ipi/2ioo) and 
V'i/2(— oo) . For pulses with Q = I the orientation of 
electron spins does not occur since such pulses rotate the 
in-plane spin components but do not generate the spin 
coherence. Quite often the pulse repetition period used 
in experiments is much shorter than spin relaxation time: 
Tr < Ts_e— • This allows to rewrite Eq. as 



2 1 + 2^i^/[r,.e(l-g2)] 



(32) 



One can see that even in the case of weak excitation 
where 1 — <^ 1 the electron spin satisfying the PSC 
reaches its utmost alignment 5.7 « — 1/2 if 



2^«1 



sin^(9/2) 
cosh^ (Try) 



(33) 



The latter equality is valid for Rosen&Zener pulses. For 
the large ratio Tg^e/Tn even quite detuned pulses, e.g., 
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FIG. 4: Electron spin polarization, , created by a train 
of Rozen&Zener pump pulses with the repetition frequency 
Tji = 13.2 ns at the moment of pulse arrival as a function 
of the electron spin precession frequency Q . Dependences 
shown on panels (a) and (b) are calculated, respectively, for 
zero detuning between the pump frequency and quantum dot 
frequency, co-p = loq , and for the detuning ujp—ujo = 1.6tt/tp . 
The chosen pulse areas are Q = 0.57r,7r and 1.5n . For g- 
factor from Ref.— the value of Q.Tu/2-k — 105 is reached at a 
magnetic field B = 1 T. The calculations used rs,e = 3 /is 
from Ref.—. 



with (wp — wo)Tp 3 , are still quite efficient in the spin 
alignment. A train of weak pulses, however, synchronizes 
electron spin precession only in a very narrow frequency 
range around the PSC. As a result the frequency depen- 
dence of S~{VL) created by such a train has a periodic 
form with sharp minima at the frequencies satisfying the 
PSC 



s: 



{nTR - 2ttNY + 

T 

$2 + (i_Q)2 + 2^(i-g) 



(34) 



Here it is assumed that 1 — Q ^ 1 , Tr ^ Tg^e and 
$ <C 1 . The width of the minima is proportional to 



For the tt -pulse excitation ( Q — > ) which can be 
realized only in the absence of detuning ( <I> = ) we 



arrive a' 
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1 cos QTr 

2 2 - cos ^ITr 



(36) 



The polarization S~ reaches its minimum value —1/2 
when the electron spin precession frequency satisfies the 
PSC condition. 

For a given pulse shape the phase <& and the factor Q 
arc interconnected and, in general, a change of the pulse 
area results in changes of both Q and <I> . Figure H] shows 
the z component of electron spin polarization, S~ , at 
the moment of pulse arrival calculated as a function of the 
Larmor frequency for different pulse areas under resonant 
excitation and for nonzero detuning. One can see from 
Fig. Ufa) that for zero detuning ( $ = ) and small pulse 
areas, S~ exhibits a sharp minimum as a function of 
Larmor precession frequency in agreement with the above 
analytical considerations. The increase of the pulse area 
transforms this minimum to the cosine-like curve. 

Figure [DJb) shows S^{il) created by an appreciably 
detuned pulse train, Wp — cjq = 1.6tt/tp . In agreement 
with Fig.[2]the effective pump power in this case is small, 
i.e., Q is close to 1 . An increase of the pump area 
from 7r/2 to tt leads to a small decrease of Q and a 
nonzero value of $ . As a result the minimum at the 
spin precession frequency satisfying the PSC becomes 
deeper and wider. The further increase of the pump area 
results in additional widening of the minimum due to the 
increase of $ but the its depth becomes smaller since Q 
starts to increase. 



III. PROBING SPIN DYNAMICS IN 
QUANTUM DOTS 

The detection of the QD spin polarization in pump- 
probe Faraday and Kerr rotation experiments is carried 
out by a linearly polarized probe pulse. The electric field 
of the probe pulse shown in Fig. 1 oscillates along the x 
axis and similarly to Eq. ([7]) it can be written as 



^v(f:) +2(l-o):^ + (l-Q)'+*^ 

i.e., it is determined either by the spin relaxation rate 
T^^^ or by the effective pump area 1 — Q and phase $ . 
Figure m shows just one period of this dependence. 

The modulation of the electron spin polarization 
S~{Cl) becomes weaker with the increasing detuning. 
For example, assuming $ = ±7r/2 and using the con- 
dition Tr <C Ts^e one can derive for arbitrary Q 



1 _^ cos ^j^^ 



2 + 2-cosnTR 
which is a much smoother function of Q . 



(35) 



£;P'-(r-,i)=i?r(r-,t)- 



V2 



c.c. 



(37) 



Here we assume that E^^{r,t) oc e~'"p'* , where 
the carrying frequency of the probe light. In order to cal- 
culate Faraday and Kerr rotation angles of polarization 
plane of the probe pulse as well as the cUipticity we first 
find the QD polarization induced by the probe field and 
then calculate the secondary electric field induced by the 
QD array. 



A. Probe-induced polarization of QDs 

Before the probe pulse arrival the electronic state of a 
QD is described by the wavcfunction We consider 
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the general case where the QD is characterized by the 
population of the electron ng = |'/'i/2p + |V'-i/2p and 
trion rit = |'03/2p + |^-3/2p states and the spin polariza- 
tion of these states Sz — (|V'i/2p ~ IV'-i/2p) /2 7^ and 
Jz = (IV'3/2p- IV'-3/2p)/2 7^ 0, respectively. Solving 
Eqs. llOl), dm) in the lowest order in we find the 

probe-induced corrections to the electron and trion com- 
ponents of the wave function 

Si^±3/2 = ^'±1/2/ ^e-'-«(*-*')di' , 

Si^±i/2 = V'±3/2 1_* ^^e'-°(*-*')dt' , (38) 

where 

Vit)^-^ldir)Erir,t)d'r. (39) 

The electron-trion superposition excited in a QD by 
the light pulses creates a local polarization, whose mag- 
nitude depends on the components of the wave function 



described in Eq. ([8]). According to the selection rules, the 
circular components of the QD dielectric polarization can 
be written as 

Pa+{r) = d*(r-)(-03/2+'^'03/2)('/'l/2+<5V'l/2)*+C.C. , 

P,-{r) = d*(r)(V'_3/2 + (5V'_3/2)(^-i/2+'5V'-i/2)* 

+ c.c. , (40) 

where the effective transition dipolc is defined by Eq. 
(fT2|l . In Eqs. (j40ll the zero-order contributions, which 
are proportional to V'3/2^1/2 ^^'^ ''/'-3/2'0!li/2 1 deter- 
mine the QD emission due to the presence of photoex- 
cited trions. They make no contribution to the measured 
pump-probe rotation signal and will not be considered 
further. The other contributions in Eqs. (|40p are induced 
by the probe pulse. The Faraday and Kerr rotation of the 
probe light polarization as well as its ellipticity are deter- 
mined only by the terms linear in Sip±i/2 and Sip±3/2 ■ 
Combining Eqs. (|40[) we can write the linearly polarized 
components of the QD polarization induced by the probe 
pulse as follows 



rie - ntr 

m 

Sz Jz 



d*(r) / dV / dtV'"°(*'-*)d(r')£;r(r',i')+c.c, 



(41) 



—00 
t 



The light wavelength is usually much larger than the size of self-organized QDs. This allows one to extract the probe 
electric field E^^{r' ,t') from the integral and present the QD polarization in the approximate 5 -function-like form 

5PS''{r.Rj,t)^5{r-R,)\lo.{R,,t) (a = a;,y), 

where Rj is the position of j -th QD. The resulting probe-field induced polarization of a single QD can be expressed 
as 



n.(H„i) = -!li-^|p|2y ^e'-(*'-*)£;r(i?,,t')dt' + c.c., 

via the integral QD transition dipole 2? = / d'^rd(r) related to the two-particle wave function F(r,r) 



(42) 
(43) 



j F(r,r)d^r 



As one can see from Eqs. (|4T|) the probe-induced QD 
polarization SP'^^ consists of two components. The 
first, 5P^^ , is parallel to the probe polarization plane 
and its magnitude is proportional to the difference of 



electron and trion occupation numbers. The second com- 
ponent, SP^^ , is orthogonal to the probe polarization 
plane and its magnitude is proportional to the difference 
of the electron and trion spin polarizations. The latter 
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polarization component is responsible for the probe-pulse 
polarization plane rotation, i.e., spin Faraday and Kerr 
effects, and for circular dichroism (ellipticity of the trans- 
mitted or reflected probe beam). Note that the appear- 
ance of the 6P^^ component is not a direct consequence 
of the spin-orbit interaction: this component is not rela- 
tivistically small as compared with SP^^ . 



Eq. (|47p can further be transformed into an integral equa- 
tion 

(t") V -r')il + q-^ grad div)Ptot (r', t) . 

(49) 



B. Circular birefringence and dichroism induced 
by photoexcited QDs 

Once the probe-induced dielectric polarization of the 
QD is known, it is possible to calculate an electric field in- 
duced by the QD ensemble and, therefore, find the probe 
polarization plane rotation and ellipticity. First we con- 
sider an experimental situation where pumping and prob- 
ing are carried out on a planar array of QDs. Then we 
generalize the results to a stack of such QD planes and a 
bulk array of QDs. 

Let us consider a layer of self-organized QDs forming 
the plane z ~ . The total electric field E— in the sys- 
tem can be represented as a sum of the incident electric 
field £?Q'^(i)e"'^ and the electric field induced by the QD 
dielectric polarization SP'^^ . The field E satisfies the 
electromagnetic wave equations 

AE{r,t)~ gT&ddivE{r,t) = -(^^)^JD, (44) 

divD = 0, (45) 

with the material equation 

Z)(r, t) = ebE{r, t) + 47rPtot(r-, t) . (46) 

Here £b is the dielectric constant of the cap layer as- 
sumed to coincide with the background dielectric con- 
stant of the QDs; Ptot{r,t) = Y.^SP^°{r,Rj,t) is the 
sum of the probe-induced polarizations over all QDs; 
D is the displacement field; ujpr is the carrying fre- 
quency; and c is the speed of light in vacuum. Al- 
though the pumping and probing of QDs is performed 
by short pulses, their duration Tp is assumed to exceed 
by far the period of electro-magnetic field oscillations 
2'K/ujpi- . Therefore, the solutions of Eqs. (|44|) - ((46l) are 
quasi-monochromatic waves with slowly varying ampli- 
tudes. 

It follows from Eq. (gS]) that div£; 
— (47r/eb) divPtot(''') which allows us to rewrite Eq. ([44]) 
in the form 

AE{r,t) + q^E{r,t) = 

_ 47r (^) ' (1 + grad div)P,,,{r, t) , (47) 

where q = t^piy^/c. By introducing the Green's func- 
tion for the three-dimensional space 
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FIG. 5: Schematic image of the light propagation in a QD 
structure, which contains a QD layer overgrown by a cap layer 
of the thickness L . Long arrows show the transmission and 
reflection of the incident probe light (with the electric-field 
amplitude Eq ) on the external surface while short arrows il- 
lustrate the creation of the secondary field, SE , due to scat- 
tering of the transmitted light by QDs. rio , tio and toi are 
the corresponding reflection and transmission coefficients. 

A plane wave of the probe electromagnetic field prop- 
agating along the z direction creates a dielectric polar- 
ization in QDs randomly distributed in the plane z = . 
Assuming the interdot distances to be smaller than the 
light wavelength one can neglect the q~^ grad div term 
in the integral of Eq. (|49p and average Ptot('»') over the 
distribution of QDs. As a result we can replace Il{Rj,t) 
by the coordinate independent vector Tl{t) . This allows 
us to rewrite Eq. in the following form 



47r 



E{r,t) ^ E^'{t)e''^''~ 

U{t)J2 I d'p'dz'G{r~r')S{p,~p')S{z'), 



(50) 

where pj are the QD positions in the two-dimensional 
layer. Assuming that QDs in the array are identical and 
randomly distributed we replace the sum in Eq. (|50|) by 
the integral and arrive at 



E 



(51) 

where Nq^ is the two-dimensional density of QDs. De- 
riving Eq. (|5ip we used the following property of two- 
dimensional integral 



(48) 



exp{iq^y z^ + p-) 1 
47r^z2 + p2 2q 
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which can be proven by adding a small positive imaginary This equation allows one to find magnitudes of the Fara- 



part to q and setting it to +0 . 

It is convenient to represent the electric field in a QD 
sample as 



day rotation signals and ellipticity. 



E{r,t) = El'{t)e"i'' ^5E{t)e 



9kl 



(52) 



If the sample contains M layers of QDs and the stack 
thickness d is smaller than the light wavelength then 
where the first term is just the incident probe field and the second term in Eq. ([52]) should merely be multiplied 
the second term describes the secondary field induced by by M . In a more conventional description of threc- 



QDs: 



dimensional (3D) ensemble of QDs, the factor MN^J^ 
can be rewritten as Nq^ , where Nq^ is the 3D con- 
(53) centration of QDs. 



The straightforward calculation shows that the Faraday rotation signal defined by Eq. ([T]) can be presented as 

/oc 
Re[E^^:it)6Ey{t)]dt, (54) 
-oo 

where only a contribution linear in SE is taken into account. Substituting SEy from Eqs. (|43|) and (|53|) we arrive at 



2 "*Q-D 



dt 

— OO ^ — oo 



di'e'-o(*'-*)£;P>;;(t)ijp;(t') 



where tqd is the radiative lifetime of an electron- hole pair confined in a QD: 

TQO 3 Ehh 

Similarly, the calculation of the ellipticity defined by Eq. ([4]) results in 

\^{El]:{t)8Ey{t)\dt . 



£ 

Substituting 5Ey from Eqs. (|43l) and (|53|) we obtain 



E 



dt 



dt'e'-(*'-*)i?p:;(ix;(t') 



(55) 



(56) 



(57) 



(58) 



In samples with a cap layer, see Fig. [SJ the Faraday 
rotation and ellipticity signals acquire an extra factor 
iioioi in Eq. (j55p . where toi and tig are the transmis- 
sion coefficients through the interface from the cap layer 
to vacuum and vice versa, respectively. 

The cap layer strongly influences the Kerr effect, i.e., 
the polarization plane rotation in the reflection geometry. 
This happens because its magnitude is determined by the 
interference between the probe beam reflected from the 
cap layer and the secondary wave induced by the QDs, 
see Fig. [5] The phase difference of the reflected and sec- 
ondary waves is determined by the cap layer thickness L 
leading to the following expression for the Kerr rotation 
magnitude 

K. = roitoi^io [cos {2qL)T + sin {2qL)E] , (59) 
where T and 8 are given by Eqs. (|55p and ([55]) . re- 



spectively, and roi is the reflection coefficient from the 
vacuum - cap layer interface. It is seen that the Kerr ef- 
fect measures, in general, a superposition of the Faraday 
rotation and ellipticity signals. 

Equations ((55| and ((58|) demonstrate that Faraday, 
ellipticity, and therefore Kerr signals (see Eg. ([55]) ) are 
proportional to the difference of electron and trion spin 
polarization in QDs: Sz ~ Jz ■ The magnitudes of the 
effects are proportional to the QD density and increase 
with a decrease of the radiative lifetime tq^ due to an 
increase of the transition dipole moment. 

In order to analyze the dependence of the Faraday 
and ellipticity signals on the detuning between probe 
frequency, , and trion resonance frequency, loq , we 
represent probe field as E'^^{t) = _E(°)s(t)e~'"p'* , where 
s{t) is the envelope function. It can be seen from 
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Eqs. ([551) and ([581) that 

T oc Im G(a;pr — wq) , £ cx Re G(ti'pi. — cjq) , (60) 
where 



G(A) = 



dt 



— CX3 ^ — CX3 



dt's(t)s(i')e 



/X iA(t-t') 



with A = Wpi- — Wo ■ It can be recast as a half axis Fourier 
transform of the probe autocorrelation function 



G(A) = 



die' 



At 



(61) and calculated for particular pulse shapes as follows 
I 



( 1 

— (iArp + 1 - e'^^p) for s{t) = 1 if - Tp/2 < t < Tj,/2 and s{t) = otherwise, 
G(A) = J ^C(2, i - for s{t) = cosh-i(7rt/rp) , 



(62) 



^2 + iATp(3 + AVp^) 



for s(t) — e^ 



where ^(a, b) is the generalized Riemann ( -function define d as ^(a, b) = X]fcln(^ + ^) ° ■ 
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FIG. 6: The dependence of real and imaginary parts of G 
on detuning, A = Wpr — wo . Dash-dot and solid curves 
are calculated, respectively, for (a) the rectangular and (b) 
Rosen&Zener pulses. 



Figure [6] shows the real and imaginary parts of G cal- 
culated for the pulses of rectangular shape (panel (a)) 
and Rosen&Zener shape (panel (b)), respectively. For 
the both pulse shape the dependences G(A) look quali- 
tatively very similar. One can see that the Faraday and 
ellipticity signals are, respectively, odd and even func- 
tions of the detuning. Therefore, the ellipticity reaches its 
maximum sensitivity for zero detuning, A = , whereas 
the Faraday signal is most sensitive for the detuned probe 
pulses with \A\Tp w 1 . The Kerr signal dependence on 
detuning has, in general, an asymmetric profile due to a 
combined contribution from and £ to its magnitude, 
see Eq. (|59l) . As shown below, this different frequency 
sensitivity leads to the different time-resolved Faraday, 
Kerr and ellipticity signals in a QD ensemble with inho- 
mogcneously broadened resonant transition energies. 



C. Effective media approximation 

At the end of this section we establish the link be- 
tween the microscopic approach developed above and 
semi-phenomenological effective medium approximation 
which is a standard tool for the description of Faraday, 
Kerr and ellipticity effects in bulk systems. We demon- 
strate below that the effective medium approximation 
can describe the Faraday and Kerr effects in a 3D en- 
semble of QDs provided that the QD density satisfies 
certain conditions. 

The irradiation of bulk homogeneous semiconductors 
with circularly polarized light creates non-equilibrium 
population of electrons and holes as well as non- 
equilibrium orientation of their spins and. thus, a 
nonequilibrium macroscopic magnetization. Hence, the 
semiconductor after the absorption of a circularly po- 
larized pump pulse changes its dielectric and magnetic 
properties. These modifications can be tested by the lin- 
early polarized probe light: the polarization plane rotates 
after its transmission through or reflection from the pho- 
toexcitcd medium leading to the Faraday and Kerr ef- 
fects, respectively. The effects are phenomenologically 
described by the Fourier component of the displacement 
field, D{uj) , which is connected with the Fourier com- 
ponent of the electric field of the probe light, E{uj) , by 



D = EbE + SiE + i[E x g] 



(63) 



Here is the background dielectric constant. Si is the 
spin-independent modification of the dielectric tensor s 
due to the filling of the conduction- and valence-band 
edge states by the photoexcitcd carriers, and g is the 
gyration vector pointing in the direction determined by 
the spin orientation S of photoexcited carriers and the 
point-group symmetry of the system. In bulk cubic semi- 
conductors g on S and the tensor Si reduces to a scalar 
Se. 
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The same description can be used for the 3D ensemble 
of QDs if their concentration Nq^ satisfies two condi- 
tions. Firstly, Nq^ should be sufficiently small so that 
the QDs may be considered as independent dipoles. Sec- 
ondly, Nq^ should be sufficiently large to have the typ- 
ical distances between QDs smaller than the light wave- 
length. The satisfaction of these conditions allows one to 
neglect the non-locality of the QD response and represent 
the displacement field D in the QD sample as 



D{r, t) = £bE{r, t) + ATT5P{r, t) , 



(64) 



where the optically induced dielectric polarization is re- 
lated to the electric field by 



6P{r,t)^ f k{t-t')E{rXW 

J — OO 



(65) 



Using Eqs. (j42|) . (j43|l we can present nonzero components 
of the tensor x(t) as 



ntr 



t3D 



\V\'e- 



(66) 
(67) 



It follows then that the QD contribution to the frequency- 
dependent dielectric permittivity tensor can be written 



h(LUo — LO — iO) ' 

-Exyiio) = (5, - J.) ^, , Q^;' .(68) 



h{ujo — cij — iO) 



Comparison with Eq. (|63[) shows that the gyration vector 
g in the photoexcited QD medium has only one nonzero 
component, 



9z 



h{ujQ — cij — iO) ' 



(69) 



which is proportional to the difference of spin densities 
of electrons and trions in the system and has a resonance 
at the trion excitation frequency. On the other hand, 
the modification of diagonal components of the dielectric 
tensor e^x = £yy is proportional to the difference in 
population of the electron and trion levels irrespective to 
the their spin orientation. 



IV. TIME-DEPENDENT TRACES OF 
PUMP-PROBE KERR AND FARADAY 
ROTATION AND ELLIPTICITY SIGNALS 

In this section we apply the derived general expres- 
sions to calculate the typical time-dependent traces of 
two color pump-probe Faraday and Kerr rotation (FR 



and KR) signals as well as the cllipticity created by short 
pulses of the resonant light and by a train of such pulses 
in an ensemble of singly charged QDs. The real QD 
structures possess two important properties that affect 
strongly the time dependent traces, but have not been 
considered in the previous sections. They are (i) inhomo- 
geneity of a QD ensemble expressed in dispersion of the 
QD resonant transition energies and electron g -factors, 
and (ii) dispersion of electron precession frequencies con- 
nected with fluctuations of the nuclear contribution to 
these frequency. Here we conduct calculations for the 
QD ensemble assuming that the broadening of trion res- 
onance frequency, dispersion of electron g -factors and 
fluctuations of the nuclear contributions to the electron 
spin precession frequency are similar to those in the sam- 
ples studied in a series of workiM^^^^^^^. Those samples 
contained 20 layers of InGaAs QDs self-organized during 
the molecular-beam epitaxy growth. 

In the following calculations wc neglect the scatter 
in the QD oscillator transition strengths and the nu- 
clear induced frequency focusing effecij^^. For illustrative 
purposes below we show the FR and ellipticity signals 
created by the electron spin polarization only. We ne- 
glect the trion -dependent contribution to these sig- 
nals in Eqs. (|55p and ([58|) which affects the time depen- 
dence traces only during the trion recombination time 
tqd ~ 400 psi^. As a result the calculated dependences 
can be directly compared with experimental data only 
for times longer than tqd . 



A. Modeling of inhomogeneities in a QD ensemble 

To model time dependences of the FR and ellipticity 
signals generated by the resonant pump pulses of circu- 
larly polarized light we assume that the distribution of 
the resonant transition energies. Popt(t^o) , hi the QD en- 
semble has the Gaussian form: 



Popt(t^o) = exp 



h^{ujQ - t^o)^ 



2(A£;)2 



(70) 



where (Dg is the average trion transition frequency and 
AE' is the half-width of this distribution. The distri- 
bution is shown in Fig. [71 In the calculations we used 
too = 1.4 eV and /^E = 6.5 mcV from Ref.^^^. Only 
a small part of this distribution is excited by the pump 
pulse with Tp = 1.5 ps. This part is proportional to the 
pulse spectral width ^ h/rp = 1.75 meV and is centered 
at the pump carrier frequency oj-p . The photoexcited 
part of the QD distribution is shown in Fig. [7] by filled 
Gaussian at low-energy part of Popt(^o) ■ 

The dispersion of the electron spin precession fre- 
quency in a QD ensemble is determined by the dispersion 
of electron g -factors and fluctuations of the nuclear con- 
tribution to the precession. The electron g -factor, g^ , 
depends generally on the effective energy gap of the QD, 
i.e. on the optical transition frequency, as well as on the 
QD shape and compositioi:^. The first effect gives rise to 
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a correlation between the average value of g -factor and 
the trion resonance frequency and can be approximated 
by a linear function 



5e(wo) = AhhJo + C . 



(71) 



where A and C are fitting parameters. This re- 
sults in the dependence of Larmor precession frequency, 
^^l('^o) = M-B5e(^o)-B/fi. , on a trion optical resonance 
frequency, which is shown in Fig. [7] by thick inclined 
line. The spread of J^l connected with the g -factor 
dependence on an excitation frequency is controlled by 
the pump pulse width h/Tp . This spread is marked by 
a green/gray segment on a linear dependence of fli^{u!o) 
in Fig. [71 The distribution of electron spin precession 
frequencies, p{i^L) created by this effect is shown in the 
inset of Fig. [7] by filled Gaussian. We use in our calcula- 
tions A = —1.75 /ieV~^ and C = 2.99 taken from the 
fit of experimental data in Ref.— . 

The frequency dependent regular part of electron g - 
factor in Ea. ((7T|) does not provide by itself a complete 
description of electron spin precession frequency disper- 
sion connected with g -factor distribution. This disper- 
sion in a QD ensemble is strongly affected by the QD 
shape and composition. The corresponding distribution 
of g -factors can be phenomenologically described by the 
Gaussian, Pg{ge) , with the root mean square of electron 
g -factor distribution. Age . 

The dispersion of electron spin precession frequencies 
is affected also by fluctuations of hyperfine fields of nu- 
clei that are collectively acting on the localized elec- 
tron in a QD. The electron spin precession frequency 
rj = + i^N = PBgeB/h + wtv contains a nuclear 
contribution, , which is proportional to the projec- 
tion of the nuclear spin polarization on the external field 
(if external field is much larger than nuclear field fluc- 
tuations, which usually is the case)^. The nuclear con- 
tribution is connected with statistical fluctuations of the 
nuclear spin polarization in a QD. The fluctuations are 
described by a Gaussian with the dispersion Awtv pro- 
portional to N~^^^ , where is the number of nuclei 
m the QD volume^. We ignore the nuclear induced fre- 
quency focusing effeclj^, which could modify the density 
of electron spin precession mode to a comb-like shape in 
a QD ensemble exposed to a pulse train excitation Ref.—. 

The resulting broadening of electron spin precession 
frequencies connected with g -factor dispersion and with 
nuclear fluctuations, p{^) — p{fiBgeB/h + lun) is also 
described by the Gaussian: 



Pin) 



27rArj 



exp 



where ni^{ujo) 



(72) 

and Ail = 

^ [pB^geB /hy + (Alun)'^ is the total frequency 
dispersion. The range of electron spin precession fre- 
quencies generated by the pulse due to this dispersion is 
shown in Fig. [7] by crosshatched region around the linear 
dependence Qi^(ujq) . The calculations were conducted 
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FIG. 7: Schematic illustration of the resonant transition en- 
ergies distribution Popt(i^o) (black solid line) and the part of 
this distribution excited by short resonant pulse (green/gray 
profile) . The range of Larmor frequencies created by the short 
pulse is shown by the crosshatched region around the linear 
dependence ^1i,{(jjo) . Inset shows the distribution of the elec- 
tron spin precession frequencies p(fi) with (blue/solid line) 
and without (filled green/gray area) shape and composition 
contribution to the g -factor dispersion and nuclear fluctua- 
tions. 



for Age = 0.0037i^ and Aujn = 0.37 GHz extracted 
from the amplitude of random nuclear fluctuation 
fleld of 7.5 mT— . We assume that Ag^ and Awjv 
is independent of the QD resonance energy. One can 
see that for used set of Age and Alon the dispersion 
of electron spin precession frequencies p{n,) is much 
broader than one created by the g -factor dependence 
on the excitation frequency (see insert in Fig. [7]). This 
additional broadening leads to the fast dephasing of 
electron spin polarization and should be taken into 
account in a description of the time dependence of FR 
and cllipticity signals. 

To obtain the time-dependent traces of the FR and cl- 
lipticity signals for the QD ensemble we average Eqs. fSS]) 
and (j58p over the distribution of optical transition en- 
ergies, /Copt(wo) , described by Eq. ((70)) and over dis- 
tribution of electron spin precession frequencies, p(f2) , 
described by Eq. (|72p . Without any calculations, how- 
ever, one can notice from Eq. (|55p that in degenerate case 



when ajpr ~ uip , FR signal vanishes if Sz and and 
all dispersion functions are even functions of the detun- 
ing, Wpr — ujo , because the signal is proportional to the 
odd function of detuning Im[G'(a;pr — wo)] ■ For degener- 
ate case FR signal could arise for excitation at one side 
of Popt(wo) distribution, or as a result of dependence of 
electron g -factors or oscillator transition strengths on 
the optical transition energy. 
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B. Effects of a single pump pulse 

We start by considering the two color FR and elliptic- 
ity signals excited by a single pump pulse as function of 
the time delay between pump and probe pulses with fre- 
quencies Wp and Wpr , respectively. To clarify qualitative 
differences between FR and ellipticity signals we assume 
here and in the subsection llV CI that Popt(wo) is indepen- 
dent of ujq or, equivalcntly, that the resonant excitation 
of QDs is performed at the maximum of Popt('^o) 7 which 
is so broad that AE ^ h/Tp . To obtain nonvanishing 
FR signal, however, we take into account the dependence 
of an electron g -factor on the resonance transition fre- 
quency, ojQ , see Eq. ([71]). 
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FIG. 8: Time resolved dependence of the pump-probe Fara- 
day rotation and ellipticity signals initiated in the QD ensem- 
ble by a single pump pulse for (a) the degenerate, Wp = Wpr , 
and (b) nondegenerate, ujp 7^ ojpr , regimes. For clarity, the 
signals are vertically shifted from each other. Calculations 
are conducted for the pump pulse with area = vr in the 
magnetic field B = 1 T by using g -factor spectral depen- 
dence described by Eq. (|7ip and the parameters Tp = 1.5 ps, 
huipT = 1.4 eV and Ts,e = 3 /is taken from Ref.— . Insets 
show the Faraday rotation and ellipticity signals in a small 
range of delay times. Faraday rotation curve in inset in panel 
(a) is multiplied by a factor of 2. Panels (c) and (d) show 
Sf and Im[G(a;pr — 1^0)] as a function of uipr — loo for the 
degenerate and nondegenerate regimes, respectively. 

Figure EJa) shows the traces of FR and ellipticity sig- 
nals for degenerate case ( (jJp = LOp^ ) calculated in a mag- 
netic field B = 1 T for the QD ensemble with the aver- 
age g -factor, its dependence on uq and Tg^e extracted 
from the data of Refs.— i^^. The trace of the ellipticity 
signal demonstrates damped oscillations with the decay 
determined by the dispersion of electron spin precession 
frequencies. 

Figure m a) shows also that the FR signal is absent at 
zero delay time, as it was expected, due to symmetric dis- 
tribution Sz and Popt(i^o) around pumping frequency. 
Surprisingly, however, this signal is growing in time. This 



happens because the electron spin distribution created 
by the pump pulse being initially symmetric around Wp 
(see Fig. ISlJc)) is gradually loosing its symmetry due to 
different Larmor precession frequencies on the low- and 
high-energy wings of the QD distribution as described by 
Eq. (j7ip . This imbalance of the electron spin polarization 
connected with electron g -factor dependence on ujq re- 
sults in the growth of the FR signal with time. The inset 
in Fig.tSja) shows also a phase shift between the Faraday 
rotation and ellipticity signals. Calculations show (not 
presented) that oscillation frequencies of the FR and el- 
lipticity signals are also slightly different. The effect is 
connected with different spectral contributions to the FR 
and ellipticity and results in a weak time dependence of 
the phase shift. 

Figure [BJb) shows FR and ellipticity signals for the 
nondegenerate case where pump and probe pulses are 
detuned: Wp — tUp^ = 0.87r/rp . The electron spin polar- 
ization created in this case is not a symmetric function of 
iVpr — uq as one can see in Fig. Hl^d), and the probe light 
measures the electron spin polarization only at one of 
the spin distribution wings. Therefore, at i = the FR 
signal is nonzero and its time dynamics is qualitatively 
similar to that of the ellipticity signal. It is seen from 
inset in Fig. [5Jb) that the phase shift between Faraday 
rotation and ellipticity signals is close to tt . The phase 
shift and the sign of the FR signal, correspondingly, de- 
pends on the sign of the pump-probe detuning because 
Im[(J'(cjpi- — ujq)] is an odd function of the detuning. 



C. Effects of an infinite train of pump pulses 

Figure [9] shows the time dependent traces of the FR 
and ellipticity signals initiated by a train of short pulses 
of circularly polarized light with the repetition period 
Tn = 13.2 ns in degenerate regime uj-p = LUp^- . Calcula- 
tions were conducted for the pulse duration Tp = 1.5 ps 
and magnetic fields B = 1 T and B = 5 T (panels (a) 
and (b), correspondingly). Panel (c) shows the results 
calculated for Tp = 100 fs and B = 1 T. The com- 
monly used repetition period of the mode-locked lasers 
Tr is about 10 ns, which is much shorter than the typi- 
cal electron spin relaxation time in a QD. As a result the 
infinite train of such pulses creates a stationary distribu- 
tion of rotating spin polarization Sz{^',t) , which modi- 
fies strongly the FR and ellipticity signals from those cre- 
ated by a single pump pulse. The traces in Figs.[Ha), (b) 
and (c) are calculated by using the steady-state values of 
the electron spin polarization defined by Eq. (j^^ . Here, 
like in the previous subsection, the dispersion of electron 
spin precession frequency is described by the function 
gei^^o) given by Eq. ((7T|) . Figures[9{d), (e) and (f) show 
the electron spin distributions at the moment right af- 
ter the pump pulse arrival. One can see that the latter 
distributions are very different from those created by a 
single pulse and shown in Figs.[8l[c) and[8l[d). 

In the case of excitation of QD ensemble by an infinite 
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FIG. 9: Time resolved dependence of the single color pump- 
probe Faraday rotation and ellipticity signals initiated in the 
QD ensemble by a train of pump pulses with the repetition 
period Tr = 13.2 ns and pulse duration Tp = 1.5 ps ((a) 
and (b)) and Tp — 100 fs (c), in the magnetic fields B = 1 T 
(a) and (c), and B — 5 T (b). Calculations are conducted 
for pump pulses with the area Q = n using g -factor spectral 
dependence described by Eq. (|7H) and the parameters hujpr — 
1.4 eV and Ts,e = 3 ^s taken from Ref.— . Panels (d), (e) 
and (f) show the distributions S'j'(a;pr — oio) and Im[G(a'pr — 
ujo)] created by the pico- and femtosecond pulse trains at 
B = 1 T and B = 5 T, respectively. 



pump pulse train a steady state distribution of electron 
spin polarization is formed in accordance with Eq. (|29p . 
The modes satisfying the PSC: n = 2'kN/Tr ( N is 
a large integer. w 100 for i? = 1 T) provide an 
enhanced contribution to the electron spin polarization. 
The sum of these mode contributions to the spin po- 
larization results in the constructive interference around 
the pulse arrival time due to the commensurability of 
the spin precession frequencies with the cyclic repetition 
frequency of the train, 2t: /Tp^ . 

The shape of the steady state distribution of electron 
spin polarization depends strongly on the number of pre- 
cession modes, which satisfy PSC and, therefore, on pulse 
duration and on a magnetic field. At the relatively weak 
field, _B = 1 T, and pulse duration Tp ~ 1.5 ps the dis- 
tribution is asymmetric because only few modes satisfy 
PSC (see Fig.[9jd)). In this case, the ellipticity and Fara- 
day rotation signals are similar to each other analogously 
to the situation of detuned pump and probe, Fig.[5I^b). 

There are more modes satisfying PSC with increase 



of a magnetic field. This is because the dispersion of 
electron spin precession frequency increases linearly with 
magnetic field but the distance between the PSC modes 
27r/T/f does not change. The density of the mode satisfy- 
ing the PSC increases and the steady state distribution of 
electron spin precession frequency at the moment of pulse 
arrival becomes more dense and symmetric (see Fig. IHUe) , 
where the electron spin polarization was calculated for 
B = 5 T). In this case the Faraday rotation signal and 
the ellipticity signal become phase shifted relative to each 
other. More importantly, the Faraday signal vanishes at 
zero delay between pump and probe pulses (Fig. EJb)) 
similar to that for the case the single pulse excitation 
with degenerate pump and probe pulses (c.f., Fig. [SJ^a)). 
The shortening of the pulse duration ( Tp = 100 fs) also 
leads to an increase of the number of modes satisfying 
PSC as it is clearly seen in Fig.[8Kf), because of the spec- 
tral width of the laser pulse increases with shortening 
of the pulse duration. The increase number of modes 
results in faster decay of ellipticity and FR rotation sig- 
nals and in vanishing of the Faraday rotation signal at 
the moment of pulse arrival. In two-color experiments 
(not shown) the ellipticity and the Faraday rotation sig- 
nal time dependent traces become similar and the phase 
shift vanishes. 



D. Effect of the pump and probe spectral position 
and the electron g -factor dispersion 

Now we turn to the effects of the spectral distribu- 
tion of the QD transition energies and of the distribu- 
tion of electron spin precession frequencies described by 
Eqs. ([701) ^-^^d G21) on the time traces of the FR and ellip- 
ticity signals created by the train of pulses with repeti- 
tion period Tr = 13.2 ns. Panels (a)-(c) in Fig. [TOl show 
FR and ellipticity time-resolved dependences for differ- 
ent spectral positions of pump and probe pulses which 
are shown on the corresponding right-hand side panels. 
The panels (a) and (b) show the traces for the same pump 
and probe carrying frequencies (single color or degener- 
ate pump-probe setup). This frequencies are tuned to the 
peak of QD distribution in Fig.fTOFa) and to its left wing 
in Fig. [TUf b). The panel (c) shows the traces of the FR 
and ellipticity signals for the case when the pump and 
probe pulse frequencies are in the vicinity of the max- 
imum of Popt(wo) but they are slightly detuned with 
respect to each other. The distribution of optical tran- 
sition frequencies in the QD ensemble, Popt(^o) i shown 
in Fig. [TO] by black curves is described by Eq. ([70]) . The 
calculation shows that the inclusion of additional disper- 
sion of electron spin precession frequencies described by 
Eq. ([7^ leads to the faster decay of both the FR and el- 
lipticity signals due to a faster dephasing of electron spin 
precession in the QD ensemble. It is worth noting that 
the FR signal amplitude vanishes for degenerate pump 
and probe pulses tuned to the maximum of the QD dis- 
tribution, similarly to what was experimentally observed 
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FIG. 10: Time resolved dependence of the pump-probe Fara- 
day rotation and ellipticity signals initiated in the QD ensem- 
ble by a train of pump pulses for the degenerate, ujp = ujpr , 
(panels (a) and (b)) and nondegenerate, lo^ 7^ ujpr , (panel 
(c)) regimes. In the right-hand side panels we show the den- 
sity of the optical transition energies and the frequency po- 
sition of the pump and probe pulses. Calculations are con- 
ducted for pump pulses of the area B = tt , the magnetic field 
B = 1 T, Tp = 1.5 ps, and Ts,e = 3 fis . 



in Rcf.— . The nonzero FR signal at zero delay time seen 
in Figs. fTOT b) and [TOT c) is related to the asymmetry in 
the QD density distribution which is revealed in the case 
of pump pulse detuned from the luq ( hCuQ = 1.4 eV). In 
the case of two color experiments, traces of the pump- 
probe FR and ellipticity signals are very similar as it is 
seen in Fig. [TUTc). 

Figure [TT] compares the Faraday rotation and elliptic- 
ity signals for laser trains with different pulse duration: 
1.5 ps (a) and 100 fs (b). Calculations that take into ac- 
count spectral distribution of optical transition energies 
in the QD ensemble Popt(<j-'o) and total Larmor frequency 
dispersion where conducted in degenerate regime for res- 
onant excitation of QDs at the maximum of Popt(wo) 
( Wp = = loq). The case of ps-excitation differs from 
one demonstrated in Fig. llOf a) only by larger disper- 
sion of electron spin precession modes Af2 = 0.9 GHz 
(time dependent traces in Fig. [TUTa) were calculated for 
An = 0.5 GHz) and, therefore, they show faster decay 
of the signals. The spectral width of pulses in the fs- 
pulse train is larger than the distribution Popt(^^o) used 
in this calculation. This effectively decreases the number 
of electron spin precession modes contributing to the FR 
and ellipticity signals explaining dcphasing decay, which 
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FIG. 11: Time resolved dependences of the degenerated 
pump-probe Faraday rotation and ellipticity signals in inho- 
mogeneous ensemble of QDs excited by trains of pulses with 
duration 1.5 ps (panel (a)) and 100 fs (panel (b)), correspond- 
ingly. Calculations are conducted for the set of parameters the 
same as in Figs.[9ja) and|9jc). Note, that the time scale used 
here is different from the time scale used in Fig. [O] 



is slightly weaker than the decay shown in Figure [Sllc). 

Surprisingly, Fig. [TT] shows that signals created by ps- 
and fs-pulse trains are very similar and the signal decays 
are almost the same. This could occur only if numbers 
of electron precession modes satisfying the PSC for both 
excitations are comparable. In the case of ps-train exci- 
tation the dispersion of electron spin precession modes is 
controlled by Ag and Alon connected with the shape 
and composition fluctuations of QDs and the nuclear field 
fluctuations, correspondingly. Due to the small spectral 
width of the ps-excitation only a small part of the elec- 
tron spin precession mode dispersion is determined by 
the frequency dependence of the electron g -factor. This 
is not the case, however, for fs-pulse train excitation, 
where the dispersion of electron spin precession modes 
have significant contribution connected with frequency 
dependence of an average g -factor g{ujQ) on the spec- 
tral width of the fs-pulses. 

Finally, Fig. [T^ shows traces of the KR signal calcu- 
lated by using Eq. ([55]l for three different thicknesses of 
the cap layer and the same parameters as in Fig. [10] One 
can see that, for particular cap layer thicknesses, the KR 
signal looks like the either FR or ellipticity signals. How- 
ever, in general the trace of the KR signal is a linear com- 
bination the FR and ellipticity signals with their partial 
contributions depending on the cap layer thickness. 

Note that L = X/2ni, (top curves in Fig. [T^ corre- 
sponds to the real thickness L = 115 nm of cap layer 
of the QD sample which was investigated in RefJ^. For 
this cap layer time dependent trace of the KR signal is 
similar to that of the ellipticity signal. 
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FIG. 12: Time resolved dependence of pump-probe Kerr rotation signals initiated in the QD ensemble by a train of pump 
pulses for the degenerate (panel (a)) and nondegenerate (panel (b)) regimes. The traces are calculated for three thicknesses of 



the cap layer 2L = X/nt, CSA/n;,, and 0.75A/n(, , where 7it 
same as in Figs. llOf a) and llUf c) 



£(, . Calculations are conducted for the set of parameters the 



V. SUMMARY 

The formalism presented here provides a complete the- 
oretical description of single- and two-color pump-probe 
Faraday or Kerr rotation and ellipticity experiments in an 
ensemble of singly charged QDs. The analytical expres- 
sions describing the electron spin polarization created by 
a circularly polarized pump pulse or by a train of such 
pulses are derived. The expressions for the magnitudes 
of the Faraday, Kerr and ellipticity signals are presented. 

The developed theory shows that the pump-probe 
Faraday rotation and ellipticity experiments measure the 
electron spin precession in slightly different subsets of 
QDs of the ensemble leading to the different oscilla- 
tion frequencies and shapes of the corresponding time- 
dependent traces. The time-dependent traces of the 
pump-probe Kerr rotation signal are linear superposi- 
tions of the Faraday rotation and ellipticity signals whose 
relative weights depend on the cap layer thickness. 

The modeling of time-dependent traces of the Fara- 
day rotation signal shows their high sensitivity to the 
inhomogeneous properties of the QD ensemble, such as 



the transition-frequency dependence of electron g -factor 
and the nuclear-induced dispersion, as well as to the ex- 
citation conditions, such as pump and probe pulse detun- 
ing, single pulse versus train of pulses excitation, and the 
pumping intensity. The pump-probe Faraday and Kerr 
rotation and ellipticity experiments can provide a com- 
plementary information about inhomogeneous properties 
of QD ensembles. 
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